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Abstract 



We present results on the production of to mesons at high transverse momenta 
in 7r~Be interactions at 515 GeV/c. The data span the kinematic range 3.5 < 
Pt < 8 GeV/c. We compare the measured cross section with expectations 
from Monte Carlo QCD generators. The relative yield of u) to vr'' mesons is 
used to extract the ratio of vector to pseudoscalar meson production for the 
first generation of quarks. 

PACS numbers: 12.38.Qk, 13.85.Ni 
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I. INTRODUCTION 



The production of mesons at large transverse momenta (px) in hadronic collisions is a 
process that can be used to study the phenomenology of Quantum Chromodynamics (QCD) 
and the parton fragmentation process. Measurements of the spectra of particles produced in 
the fragmentation of partons influence our understanding of high-p^ processes and affect the 
design and implementation of Monte Carlo event generators. Since mesons constitute the 
majority of the particles produced in hadronic interactions at high-p-p, such measurements 
are of interest across a broad spectrum of experiments in high energy physics. Although 
much is known about pion production in this regime, the hadroproduction of u mesons has 
not been as extensively studied [0-0]. A comparison of uj to vr" production, which reflects 
the overall ratio of vector meson to pseudoscalar meson production {V/P), can also be used 
to sharpen the value of the V/P phenomenological parameter used in current Monte Carlo 
event generators. 

We report on uj production in 7r~Be collisions at 515 GeV/c, as measured in E706, an 
experiment which was designed to study the production of direct photons, neutral mesons, 
and associated particles at high-p^ using the Meson West Spectrometer at Fermilab |Q. The 
apparatus included a charged particle spectrometer consisting of silicon microstrip detectors 
in the target region and multiwire proportional chambers and straw tube drift chambers 
downstream of a large aperture analysis magnet 0. The target consisted of two 0.8 mm 
thick copper foils, immediately upstream of two cylinders of beryllium, one 3.7 cm long 
and the other 1.1 cm long. Photons were detected in a 3 m diameter, lead and liquid- 
argon, sampling electromagnetic calorimeter (EMLAC), located ^ 9 m downstream of the 
target 0. The EMLAC readout was subdivided azimuthally into octants, each consisting 
of interleaved, finely segmented, radial and azimuthal views. The radial views were also 
used to form a fast high-p^- event selection trigger. Trigger decisions were based on global 
(full octant) and local (sixteen 5.5 mm strips) sums of EMLAC energy weighted to measure 
transverse momentum 0. 

II. DATA ANALYSIS 

The data sample corresponds to an integrated luminosity of 8.6 events/pb. The mea- 
surement of the UJ meson cross section was based on its decay u 7r°7. Although this mode 
has a branching ratio of only 8.5% [0, its decay chain to three photons yielded a clear sig- 
nature. In the following section, we briefly describe the analysis procedures used to extract 
the u signal in these data. A more detailed description of the procedures used to select and 
reconstruct the events, and to correct the data for losses due to trigger inefficiencies and 
selection requirements, can be found elsewhere 1^,0. 

A. Signal extraction 

The invariant mass distribution for tt^'j pairs, subject to only minimal kinematic cuts, is 
illustrated in Fig. 0. In our u study, we defined combination of two photons with 

invariant mass, M^^, in the range 110 MeV/c^ < M^^ < 165 MeV/c^ and energy asymmetry 
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[A = \Ei — E2\/{Ei + £"2)] less than 0.75. The uj signal in Fig. |l] consists of a shoulder riding 
on a steeply falling background. To improve the signal-to-background ratio, we investigated 
the cos^* distribution, where 9* is defined as the decay angle of the 7r° in the 7r''7 rest frame, 
relative to the 7r*'7 line of flight. Monte Carlo studies showed that the cos^* distribution 
for accidental pairings of vr^'s and 7's is peaked near ±1 p[, whereas the distribution for 
unpolarized w's is isotropic. (Since we detected no evidence for uj polarization in our data, 
we assumed for the purposes of the analysis presented in this paper that the signal was 
dominated by unpolarized w's.) The requirement | cos^*| < 0.6 eliminates a large fraction 
of the combinatorial background |]^, resulting in the distribution shown in the insert in 
Fig. |I], in which the uj signal stands out much more clearly. 

Figure |^ displays the weighted mass distribution of 7r°7 pairs in the vicinity of the uj 
mass. These data have also been subjected to additional kinematic constraints including a 
7r° sideband subtraction which Monte Carlo studies showed to be effective in removing 777 
combinations that peaked in the vicinity of the uj signal The effect of these and other 
kinematic criteria was estimated using the Monte Carlo simulation described below. 

Fits to the 7r°7 mass spectra for different intervals of pt were used to extract the contri- 
bution of the UJ signal as a function oipx- A Gaussian shape for the signal, combined with 
a third order polynomial for the background, yielded a reasonable description of the data. 



B. Monte Carlo 

A full event simulation was used to evaluate most corrections to the cross section. This 
simulation relied on the HERWIG |TD| event generator and a geant ITT] Monte Carlo sim- 



ulation of our apparatus |p. Events containing leading uj mesons were generated and re- 
constructed using our standard reconstruction package The ratio of reconstructed to 
generated events was parameterized as a function of px-, fitted with a Theta function con- 
volved with a Gaussian (Fig. ||), and applied to the uj spectrum. Due to the relatively large 
opening angle for photons in uj decay, and the resulting spatial spread of these photons in 
the EMLAC, the decay products infrequently satisfied the E706 high-p^ trigger. The overall 
detection efficiency for the uj (in the 7r°7 decay mode) was ~ 20% at high-p^^, and differed 
between the global and local triggers at lower-p^- 



C. Systematic uncertainties 

We estimate the combined systematic uncertainty in the uj cross section due to recon- 
struction efficiency and normalization to be 15% [If. The energy response of the electro- 
magnetic calorimeter was calibrated using vr*^, t], uj, and J /if) mass peaks. The energy scale 
uncertainty was determined to be less than 0.5% 0, and contributed between 5% and 12% 
(as a function of pt) to the systematic uncertainty of the uj cross section. 

The uncertainty in the fitted background was estimated by using different parameteri- 
zations for the background, and by varying the 7r°7 mass range and the bin sizes used in 
the fit. We estimate a 10% contribution to the systematic uncertainty in the determination 
of the signal due to the fitting procedure. The uncertainty in the trigger corrections was 
estimated by comparing the ui cross section obtained using samples selected with different 
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triggers. This uncertainty ranged from 17% for uo mesons with px of 3.5 GeV/c to 2% for 
8 GeV/c. These quoted systematic uncertainties do not exphcitly incorporate contributions 
due to the possibihty of u polarization, however, as aheady stated, we detected no evidence 
for uo polarization over the range | cos^*| < 0.6. 

The systematic uncertainties, combined in quadrature, are quoted with the cross sec- 
tions in Table |. The overall systematic uncertainty on the uj cross section was 30% at 

=3.5 GeV/c, 22% at 5 GeV/c, and 20% at 8 GeV/c. 



III. RESULTS 



Table | lists our measured inclusive invariant cross section for uj meson production in 
TT^Be interactions at 515 GeV/c along with statistical and systematic uncertainties. The 
results are binned in pt from 3.5 to 8 GeV/c, and averaged over the range of our acceptance 
in rapidity (—0.5 < ycm < 0.75). The uj cross section is also displayed in Figure H, and 



compared with expectations from PYTHIA [|12| and HERWIG [0 . The predictions from both 



Monte Carlos are substantially smaller than the measured uo cross section. Comparison 
between the relative yields of cu's originating in the Be and Cu target materials resulted in 
the value a = 1.12 ± 0.07 ± 0.07 |^ using the parameterization a a oc A". The Monte Carlo 
results have been adjusted for this nuclear effect. 

Figure ^ and Table | display the relative yields of uj and 7r° mesons measured in E706, 
in terms of the ratio of the inclusive differential cross sections as a function of px- The 



prediction from the PYTHIA generator is consistent with our measured ratio. The ratio 
from HERWIG is much smaller than both our measurement and the result from PYTHIA. 
We include for comparison three previous results on the uj to 7r° cross section ratio for 
incident protons on p. Be, and C targets These earlier measurements were integrated 

over Pt and displayed at their minimum pt value. 

The UJ to 7r° ratio can be used to determine the value of V/P. The quark content of 
both UJ and vr*^ is the same, rendering their production ratio insensitive to beam and target 
composition. Corrections to V/ P to account for indirect production were determined using 
both PYTHIA []12| and HERWIG [jl3[. In both cases, the data from E706 were found to require 



only a relatively small correction for indirect contributions from the decay of higher mass 
hadrons. We note, however, that in pythia almost all uj mesons are produced directly, 
which seems to us an extreme assumption, and the HERWIG uj to 7r° ratio differs greatly 
from our data. Nevertheless, if we use these Monte Carlos to correct our measured uj to 
vr" ratio for indirect production, the resulting V/P values are 1.2 ± 0.1 using pythia and 
0.9 ± 0.1 using HERWIG. 
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FIG. 1. Invariant mass distribution of the TT 7 system subject to minimal selection criteria. The 
insert shows the corresponding mass distribution with the additional requirement of | cos^*| < 0.6. 
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FIG. 2. Invariant mass distribution of tt 7 system with the cuts and corrections described in 
the text. The background is parameterized by a third-order polynomial. 
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FIG. 3. Efficiency for detecting the co signal in the tt 7 decay mode, as a function of 
including losses due to geometric and trigger acceptances. 
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FIG. 4. Inclusive invariant differential cross section per nucleon for to production as a function 
of pt- The curves correspond to expectations from Monte Carlo QCD generators, and have been 
adjusted for nuclear effects. 
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FIG. 5. Ratio of the invariant cross section for cu production to the tt cross section as a 
function of pT- Results of this experiment are compared with previous measurements, and with 
expectations from PYTHIA and HERWIG (assuming uj and tt^ production have the same nuclear 
dependence). Results from earlier experiments are integrated over pT, and displayed at their 
minimum pT values. 
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TABLES 



TABLE I. Invariant differential cross section per nucleon [EdPa/dp^) for the inclusive reaction 
7r~ + Be — > a; + X, averaged over the rapidity range, —0.5 < ycu < 0.75. Also included is a ratio 
of the ui cross section to the corresponding 7r° cross section measured by E706. 



Pt 




Ed^a/dp^ 




(GeV/c) 


(GeV/c) 


pb/{GeY/cf 


ratio 


3.5-4.0 


3.73 


12500±2400±3800 


0.70±0.13±0.20 


4.0 - 4.3 


4.16 


4900±1000±1300 


0.89±0.18±0.22 


4.3 - 4.7 


4.49 


1910± 390± 480 


0.90±0.18±0.20 


4.7-5.0 


4.83 


840± 190± 190 


1.05±0.24±0.21 


5.0 - 5.5 


5.23 


380± 70± 80 


1.07±0.20±0.20 


5.5 - 6.0 


5.73 


130± 30± 30 


l.lit 0.3± 0.2 


6.0 - 6.5 


6.23 


75± 16± 16 


2.0± 0.4± 0.4 


6.5 - 7.0 


6.73 


22± 9± 5 


1.6± 0.6± 0.3 


7.0 - 8.0 


7.11 


1± 3± 1 


1.2± 0.9± 0.3 
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